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The Cosmic Microwave Background radiation is a fundamental prediction of Hot Big 
Bang cosmology. The temperature of its black-body spectrum has been measured 
at the present time, Xcmbr,o — 2.726±0.010 K, and is predicted to have been higher 
in the past. At earlier time, the temperature can be measured, in principle, using 
the excitation of atomic fine structure levels by the radiation field. All previous 
measurements however give only upper limits as they assume that no other significant 
source of excitation is present. Here we report the detection of absorption from the 
first and second fine-structure levels of neutral carbon atoms in an isolated remote 
cloud at a redshift of 2.33771. In addition, the unusual detection of molecular 
hydrogen in several rotational levels and the presence of ionized carbon in its excited 
fine structure level make the absorption system unique to constrain, directly from 
observation, the different excitation processes at play. It is shown for the first time 
that the cosmic radiation was warmer in the past. We find 6.0 < Tcmbr < 14 K at 
z = 2.33771 when 9.1 K is expected in the Hot Big Bang cosmology. 

One of the firm predictions of standard Big Bang model is the existence of relic radiation from 
the hot phase the Universe has experienced at early timesS. The cosmic microwave background 
radiation (CMBR) has been discovered serendipitously by Penzias & Wilsoni in 1964. The fact 
that its spectrum follows with remarkable precision a Planckian distribution over several decades 
in frequency is a strong argument in favor of the Hot Big Bang cosmology. However the presence 
of the radiation at earlier time has never been proven directly. The temperature of its black-body 
spectrum is predicted to increase linearly with redshift Tcmbr(-z) = Tcmbr(0)x(1 + z) and its 
local value has been determined very accurately by the Cosmic Background Explorer (COBE) 
Tcmbr(O) = 2.726±0.010 iJ. Detecting the presence of relic radiation at earlier epochs and 
confirming the well defined temperature evolution is therefore a crucial test for cosmology. To 
perform this test, one can use the possibility that excited fine structure levels of the ground-state 
of atomic species can be partly populated by the Cosmic Microwave Background Radiation when 
the energy separation of the levels is similar to the energy at the peak of the radiation energy 
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distribution. 

The relative populations of excited levels can be measured from the absorption lines seen in the 
spectrum of distant quasars and used therefore to constrain Icmbr a t high redshift§@i'3i'i0. 
The expected values make neutral carbon, C°, particularly well suited for this purpose. The 
ground term is split into three levels (J = 0, 1, 2) with J = 0-1 and J = 1-2 energy separations 
kT corresponding to, respectively, T = 23.6 and 38.9 K. However the fine structure levels can also 
be excited by collisions (mostly with electrons and hydrogen) and by UV pumping and following 
cascades. The ionization potential of neutral carbon, 11 eV, is below the ionization potential of 
hydrogen, 13.6 eV, and C° can only be seen in dense, neutral and highly shielded g&M. Therefore, 
excitation by collisions cannot be neglected. 

To constrain TcmbR) the kinetic temperature, the particle density of the gas and the local UV 
radiation field must be known. As it is very difficult to disentangle these various excitation 
processes, all measurements up to now have lead to upper limits on Tcmbr only. 

Our measurement has been obtained in a unique absorption system where absorption lines of 
neutral carbon in the three fine-structure level of the ground term, C°, C *, C ** are observed 
together with absorption lines of singly ionized carbon in its excited fine-structure level, C + *, 
and absorption lines of molecular hydrogen (H2) in the J = to 5 rotational levels. The 
population and depopulation of the first excited rotational level of H2 (J = 1) from and to the 
ground state (J = 0) is controlled by thermal collisions. Therefore the excitation temperature 
T01 is approximately equal to the kinetic temperature. The fine structure upper level of the 
C + ground-state doublet is mostly populated by collisions and depopulated by radiative decay. 
Therefore, once the temperature is known, the particle density can be derived from the C + */C + 
ratio. Finally, the UV radiation flux can be constrained from the populations of the J = 4 and 5 
H2 rotational levels. 



Observations 



We have used the Ultra-violet and Visible Echelle 




- 4 - 



ESO KUEYEN 8.2 m telescope at the Paranal observatory on April 5 and 7, 2000 to obtain a 
high-spectral resolution spectrum of the z cm = 2.57 and my = 18.4 quasar PKS 1232+0815. 
Standard settings have been used in both arms of the spectrograph. Wavelength ranges were 
3290-4519 A in the blue; 4623-5594 and 5670-6652 A for the red chips. The slit width was 
1 arcsec and the CCDs were binned 2x2 resulting in a resolution, AA/A, of ~45000. The exposure 
time was 3 hours in seeing conditions better than 0.8 arcsec full width at half maximum. The 
data were reduced using the UVES pipeline in an interactive mode. The pipeline is a set of 
procedures implemented in a dedicated context of MIDAS, the ESO data reduction package. The 
main characteristics of it is to perform a precise inter-order background subtraction for science 
frames and master flat-fields, and to allow for an optimal extraction of the object signal rejecting 
cosmic ray impacts and performing sky-subtraction at the same time. The reduction is checked 
step by step. Wavelengths are corrected to vacuum-heliocentric values and individual ID spectra 
are combined together. This resulted in a S/N ratio per pixel of 10 around 3700 A and 20 around 
6000 A. Typical errors in the wavelength calibration is ~0.5 km s _1 . 

Fit to the lines, metallicity and molecular content 

This is the first time C°, C°* and C ** absorption lines are detected at high redshift (see Fig. |l|). 
The detection of these three species at z a b s = 2.33771 is confirmed by the presence of several 
transitions. C° absorption lines with rest wavelengths 1139.79, 1157.91, 1260.73, 1560.31 and 
1656.92 A , C°* absorption lines with rest wavelengths 1194.40, 1279.89, 1329.09, 1329.10, 1329.12, 
1560.68, 1560.71, 1656.27, 1657.38 A and 1657.91 and C ** lines at 1657.00 and 1658.12 A are 
clearly detected and are free from any blending with absorption due to other systems. We use 
the oscillator strengths given by Welty et al.@ and standard Voigt profile fitting to determine 
the amount of matter lying on the line of sight and characterized by column densities, N, and 
the Doppler parameters, b = J b^ h + 6^ urb where 6 t h = \J2KT jm is the thermal broadening due 
to the temperature of the gas, with T the kinetic temperature, K the Boltzman constant and m 
the mass of the particle and &turb is the turbulent broadening due to macroscopic motions of the 
gas. It can be seen in Figs. |l| and ^ that the system is dominated by a single and well defined 
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component. We consistently impose the same b value to this component for all species and obtain 
a best fit for b = 1.7±0.1 kms" 1 . This value is small compared to the spectral resolution of the 
data but is ascertained by the relative optical depths in the numerous lines with very different 
oscillator strengths. Results are given in Table 1 and the fit to a portion of the spectrum is shown 
in the bottom panel of Fig. |l| and in Fig. 

The neutral hydrogen column density in the system, log N(H°) = 20.90±0.10, has been obtained 
by fitting the damped Lyman-a line. At such a high column density, hydrogen is mostly neutral 
and the elements of interest here, iron, magnesium, silicon and carbon are mainly in the singly 
ionized state. Column densities are derived from simultaneous Voigt profile fitting of all available 
absorption lines using the same Doppler parameter for all species and the same column density 
for each species. We concentrate however on weak lines as Fe + AAA1125, 1608, 1611, Si + A1808 
or Mg + A1239 because they are optically thin. Indeed in this case the strength of the line does 
not depend on the Doppler parameter and directly gives the column density. Profiles of a 
few transitions are shown on Fig. || and results are given in Table 1. The a-chain elements 
magnesium and silicon have similar abundances, Z(K) = A r (X)/A r (H), relative to solar within the 
measurement uncertainties, log Z/Zq ~ —1.2. However we note that, relative to what is seen in 
the sun, iron is underabundant by about a factor of 5 compared to silicon and nickel is slightly 
underabundant compared to iron. If these differences are due to iron and nickel being depleted 
into dust-grains, the depletion is of the same order as what is observed in warm halo gas in our 
GalaxyB. All this suggests that the system is very much like typical damped Lyman-a systems, 
with metallicity about one tenth of solar and small depletion of iron type element into dusti-0. 

Our echelle data confirms the presence of molecular hydrogen first detected by Ge and BechtoldS. 
At the resolution of the spectrum, most of the H2 lines are free from blending. The wide 
wavelength coverage of the spectrum implies that absorption lines from various Lyman-bands are 
seen. Column densities of H2 in different rotational J levels are obtained by simultaneously fitting 
various Lyman-bands which are free from contamination due to intervening Lyman-a absorption 
from the diffuse intergalactic medium. Absorption profiles corresponding to transitions from J > 1 
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rotational levels are consistent with a single component at the same redshift z^bs = 2.33771 as 
the C° lines. Lines from the J = and 1 levels are broader and most of the absorption features 
have similar profiles suggesting the presence of several components. We obtain a best fit to these 
absorption lines with four components. Results of Voigt profile fitting are given in Table 2. 
The total H2 column density is iV(H2) = 1.52xl0 17 cm -2 and the molecular fraction is / = 
2iV(H 2 )/(2iV(H2)+iV(H )) = 3.8xl0~ 4 which is about two orders of magnitude less than what 
was derived from low dispersion data. 

Physical conditions in the gas 

In molecular clouds, the kinetic temperature of the gas is given by the excitation temperature, 
T01, measured between the J = and 1 levels. For the components at z a \> s = 2.33735, 2.33754 
and 2.33793 we find Tqi = 74 ± 7, 64±7 and 66±10 K respectively. None of these components 
show detectable absorption due to C°. The characteristics of the gas are very similar to what has 
been measured recently in diffuse gas in the halo of our Galaxy0. For the z a b s = 2.33771 cloud, 
^kin = Tqi = 185±100 K. The excitation temperature for other levels, J = 2 to 5, is close to 
< T ex > = 400 K. This means that processes other than collisions are at play in determining the 
relative populations of the different levels. Indeed, the J = 4 and 5 levels can be populated by 
cascades following UV-pumping and H2 formation. We can evaluate the UV-pumping rate from 
J = to J = 4, /?(0), using the simple model of Jura§ and assuming that molecules are formed 
on dust grains with formation rate R, 

p 4 ,o/S(0)n(H 2 , J = 0) + O.WRn(H)n = A(4 — > 2)n(H 2 , J = 4) (1) 

where n = n(H) + 2n(H2), p^o = 0.26 is the pumping efficiency from level J = to level 
J = 4. We use the spontaneous transition probabilities A(A — > 2) = 2.8 x 10~ 9 s" 1 as given by 
Dalgarno & Wright!]. Scaling the formation rate of H2 in our Galaxy (R ~ 3 x 10 s cm ) 
with the dust to hydrogen ratio measured in the system (< 10~ 13 the Galactic one) it is easily 
seen that populating the J = 4 level following formation of a molecule on dust grain is negligible 
for densities less than 5xl0 3 cm -3 . We estimate the density from the excitation of C + to be at 
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least a factor of hundred less than this value (see below). Moreover, for such large densities, the 
J = and J = 2 levels as well as the J = 1 and J = 3 levels should be in thermal equilibrium 
which is not the casefi We therefore can neglect populating of J = 4 and 5 levels by cascades 
after the formation of a molecule. We finally derive from Equation 1 that the photo- absorption 
rate in the Lyman and Werner bands in the surface of the cloud is of the order of 2xl0~ 10 s _1 



which is quite modest and similar to that observed along sight-lines in our Galaxylisl. 

The excited level of the C + ground-state term is populated by collisions with electrons 
and hydrogen atoms0. Electrons cannot be neglected because the corresponding collisional 
cross-section is large. To estimate the electron density we use the ratio N(Mg + ) / N(Mg°) = 178 
which is fairly well determined from the detection of weak lines Mg°AA1737,1827 and Mg + A1239. 
As we have shown above that the ionizing flux in the cloud is similar to that observed in our 
Galaxy we can estimate the electron density by equating the Galactic ionizing rate of Mg° to the 
recombination rate of Mg + lll. We find n c ~ 0.02 cm~ 3 . Note that the electronic density cannot be 
much smaller than this value because, from the excitation of the J = 4 and 5 H2 levels, we know 
that the ionizing flux is not smaller than the Galactic value. 

Then, the hydrogen density can be derived from the A r (C + *)/A r (C + ) ratio. Unfortunately, 
the C+A1334 absorption line is saturated and cannot be used directly to infer iV(C + ). We 
therefore consider that the carbon metallicity can be derived from the silicon metallicity with 
some correction. Indeed, it is known that the metallicity of the a-chain elements is enhanced 
compared to carbon by a factor of about twoll when the mean metallicity is low (i.e. < —1.0). 
To be conservative, we consider that Z(C) = Z(Si)/2 and that silicon is not depleted into dust. 
Note that both assumptions maximize the C + */C + ratio and, as a consequence, maximize the 
derived hydrogen density, nn- In addition, the very small C°/C + value shows that the ionization 
correction is negligible. Note also that absorptions due to Si 3+ and C 3+ are unusually weak in 
this damped system. The fit of the C + *AA1335.6,1335.7 lines is shown in Fig. |2[ As can be seen 
from this figure the observed C + * profile suggests the presence of possible extra components. Our 
measured column density is therefore an upper limit which will turn into an upper limit on the 
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hydrogen density. It must be noted that the strongest constraint comes from the 1335.6 line in the 
blue wing of the profile. 

Results on the determination of nH are presented in Figure ||. For a kinetic temperature in 
the range 85 < T <285 K, as suggested from the H2 absorption, we find 20 < ran < 35 cm -3 (solid 
curve on Fig. |3j). All assumptions above maximize the hydrogen density which can therefore be 
considered for each temperature as a conservative upper limit. 

Cosmic Microwave Background radiation temperature 

The fine structure levels of C° can be populated by several processes, mainly collisions with 
hydrogen atoms and electrons and pumping due to the local UV-radiation and to the CMBR. The 
different contributions can be estimated once the temperature and the particle density are known. 
We have derived above 85 < T < 285 K, 20 < n H < 35 cm -3 and n e /n H ~ 0.001. 

Following Meyer et al.i and using cross-sections given by Nussbaumer & Ruscai and Launay 
& Roueffi, we investigate the fine-structure excitation of C° keeping the particle density and 
the temperature within the ranges determined in the previous Sections. For a radiation field 
similar to that of the Milky- wayEl, the UV pumping rate from the ground state to the excited 
states of C° is 7.55xl0 -10 s^ 1 . The corresponding value for the hydrogen collisional rate in the 
density and temperature ranges considered here is of the order of ~ 1.2 — 1.4 x 10~ 8 s _1 . Thus 
the contribution due to UV pumping is negligible. Keenan et alS noticed that collisions with 
electrons are unimportant for kinetic temperatures in the range 100-500 K and a ratio of electron 
density to hydrogen density less than 0.01. Collisions with electrons can therefore be neglected. 

The observed 2a range for N(C°*)/N(C°) and N(C°**)/N(C°) are, respectively, 0.28-0.48 and 
0.020—0.060. Assuming collision by hydrogen atoms is the only process at work, we estimate the 
range of hydrogen density needed to explain the C° populations. This density range is shown as 
a shaded region in Fig. |3|. From this figure it is apparent that the observed upper limit on the 
hydrogen density derived above is not sufficient to populate the excited fine-structure levels of C° 
and that an additional source of excitation is needed. The only process we are left with is direct 
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pumping due to photons from the relic background radiation. 

More generally, we estimate the allowed range for the CMBR temperature as a function of kinetic 
temperature, after taking into account all previously discussed processes. The allowed area in the 
TcMBR~TKinetic plane is shown as the shaded region in Fig. |4|. This region is obtained using the 
upper limit on the hydrogen density derived above and the 2a ranges of the C°*/C° and C°**/C° 
ratios. Thus the lower Tcmbr boarder of this area gives a stringent lower limit on the CMBR 
temperature at z a b s = 2.33771. 

Note that an upper limit on Tcmbr is obtained assuming that the CMBR is the only excitation 
process at work. This limit is shown as a dotted line in Fig. ||. 

The presence of a relic radiation field at any redshift is more or less admitted by most of the 
practicing cosmologists. Its reality is proved observationally only in the local universe through 
direct measurements!!!! and the Sunayev & Zeldovich effect^. Our measurement demonstrates 
for the first time that the cosmic radiation exists at earlier times and has a higher temperature 
than today. In the standard Big-Bang cosmology the cosmic microwave background is a relic 
radiation left over from an early hot phase. In such a model the radiation evolves adiabatically in 
the expanding universe and the temperature at any redshift is simply T{z) = T{z = 0)(1 + z). We 
summarise all the available upper limits in Fig. |5[ together with our measurement. The prediction 
for the adiabatic expansion is shown as a dotted line. All the upper limits and our measurement 
are consistent with the predictions of the standard model. Doing similar analysis over a large 
redshift range will provide a direct and model independent measure of the evolution of the cosmic 
microwave background radiation. 
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TABLE 1: Heavy elements from the damped Lyman-a system 



Ion 


log N (cm 2 ) 


b (kms- 1 ) [Z/H] 


[Z/H]-[Z/H] ( 


H° 


20.90±0.10 






C° 


13.86±0.22 


1.70±0.10 




C°* 


13.43±0.07 








12.63±0.22 






C+* 


< 14.00 






Mg° 


13.19±0.09 






Mg+ 


15.44±0.09 


-5.46±0.13 


-1.04±0.13 


Si+ 


15.24±0.11 


-5.66±0.19 


-1.20±0.10 


Fe+ 


14.68±0.08 


-6.22±0.13 


-1.73±0.13 



N is the column density and b the Doppler width of the line; b = (2kT/m + b^ urb ) 1//2 , where 
m is the atomic mass, T the temperature and &turb the characteristic turbulent velocity of the 
gas. [Z/H] = log iV(Z)-log N(R i) is the metallicity and [Z/H]- [Z/H] the metallicity relative 
to the value measured in the Sun. We used the solar metallicities from Savage & SembachB. All 
wavelengths and oscillator strengths are from MortonEl and Welty et al.0. 
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TABLE 2: Fit results to rotational levels in the vibrational ground state of H 2 





Level 


log N (cm 2 ) 


b (km s l ) 


2.33735 


J = 


3.60 ±0.34 x 10 15 


24.07 




J = 1 


3.34 ± 0.40 x 10 15 


24.07±10.20 


2.33754 


J = 


4.10 ±0.48 x 10 15 


14.09 




J = 1 


2.65 ± 0.44 x 10 15 


14.09±4.57 


2.33771 


J = 


2.30 ± 1.00 x 10 16 


4.62±0.36 




J = 1 


8.30 ± 2.60 x 10 16 


4.62 




J = 2 


1.59 ±0.27 x 10 16 


4.62 




J = 3 


9.39 ± 1.10 x 10 15 


4.62 




J = 4 


4.42 ± 0.80 x 10 14 


4.62 




J = 5 


4.64 ±0.77 x 10 14 


4.62 


2.33793 


J = 


1.00 ±0.04 x 10 15 


21.74 




J = 1 


6.69 ± 0.40 x 10 15 


21.74±1.40 
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5528 5530 5532 5534 

Observed wavelength (A) 

Fig. 1. — A sample of H2 and C° absorption lines at z a b s = 2.33771: Portions of the 
normalized spectrum of the quasar PKS 1232+0815 taken with the Ultra-violet and Visible Echelle 
Spectrograph mounted on the 8.2 m KUEYEN telescope of the European Southern Observatory 
on the Paranal mountain in Chile. Upper panel: A selection of H2 absorption lines from the J = 0, 
1, 2 and 3 rotational levels from the v = 0-1 Lyman band. The model fit with parameters given 
in Table 2 is overplotted to the data as a dashed line. Lower panel: Detection of absorption lines 
from C°, C°* and C ** at z ahs = 2.33771 in the damped Lyman-a system. The model fit with 
parameters given in Table 1 is over-plotted to the data as a dashed line. 
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Relative velocity 



Fig. 2. — A sample of heavy element absorption lines at z a ^ s = 2.33771: Absorption 
profiles of a few transitions (indicated on each panel) from the damped Lyman-a system toward 
PKS 1232+0815. The normalized flux is given on a velocity scale with origin at z = 2.33771. The 
C + * profile is slightly broader than the other lines. The C + * line is a blend of C + *A1335.71 and 
C + *A1335.66. Fitted models are over-plotted as dashed lines. The fit is performed over all the 
absorption lines available in the spectrum. Our best fitted column density of 10 14 cm~ 2 for C + * 
overpredicts the absorption at A1335.6. 
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Fig. 3. — The hydrogen density as a function of kinetic temperature in the z a } 3S =2. 33771 
cloud: The continuous curve with downward arrows shows the strict upper limit on the hydrogen 
density estimated from the fine-structure excitation of C + . The shaded area gives the 2a range 
for the hydrogen density required to explain the population ratios in the fine-structure levels of 
C° assuming there is no cosmic background radiation. All other excitation processes (collisions 
by electrons and hydrogen atoms, UV pumping) have been taken into account. It is apparent 
that the observed density is too small to explain the C° excitation demonstrating that the Cosmic 
Microwave Background exists at z = 2.33771. 
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Fig. 4. — Cosmic microwave background temperature as a function of kinetic 
temperature of the gas: For a given temperature we obtain an upper limit of the hydrogen 
density from the C + fine-structure population. This is used to derive the contribution of hydrogen 
collisions to the C° fine-structure excitation. The excess is used to determine the range for the 
CMBR temperature. The shaded region gives the 2a range of the Cosmic Microwave Background 
Radiation temperature allowed by the observed population ratios of neutral carbon fine-structure 
levels. The horizontal dotted line is the upper limit on Tcmbr assuming CMBR to be the only 
excitation process. The dashed line is the predicted temperature at z = 2.33771 in the standard 
big— bang model. This demonstrates the presence of a background radiation with a temperature 
at least twice that measured in the local universe, Tcmbr o = 2.726 K. 
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Fig. 5. — Measurements of the Cosmic Microwave Background Radiation temperature 
at various redshifts: The point at z = shows the result of the Cosmic Background 
Explorer (COBE) determination!, Tcmbr(O) = 2.726±0.010 K. Upper limits are previous 
measurement using the same techniques as in this paper. We also include our two new 

unpublished upper limits at z = 2.1394 along the line of sight toward Tololo 1037—270. The 
measurement from this work, 6.0 < Tcmbr < 14.0 K at z = 2.33771, is indicated by a vertical bar. 
The dashed line is the prediction from the Hot Big Bang, Tcmbr(z) = 2cmbr(0)x(1 + z). 



